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Introduction 

Conventional  Fourier  series  analysis  enables  one  to  extract  information 
from  a periodic  function  of  time.  In  practical  situations,  however,  one 
often  deals  with  finite  time  functions,  where  conventional  Fourier  series 
techniques  are  inadequate.  Finite  functions  result  from  a singular  or  finite 
set  of  disturbances.  The  transient  response  of  a circuit  is  an  example  of 
such  a disturbance. 

Often  one  deals  with  network  synthesis  problems  where  it  is  desirable 
to  describe  the  output  and  input  functions  of  the  unknown  networks.  Laplace 
transform  techniques  often  fail  when  the  output  of  the  network  is  a finite 
function;  this  is  because  the  Laplace  transform  of  the  finite  function  results 
in  a transcendental  function  in  the  complex  frequency  variable.  Further,  one 
cannot  express  the  finite  output  as  a conventional  Fourier  series  to  obtain  a 
rational  Laplace  transform  of  the  output. 

" The  theory  incorporated  in  finite  Fourier  filters  enables  one  to  deal 
with  the  type  of  time  functions  mentioned  above.  This  technical  note  attempts 
to  summarize,  in  a logical  development,  the  basic  theory  of  finite  Fourier 
analysis,  and  the  synthesis  which  incorporates  the  theory  in  the  finite 
Fourier  filters. 


Finite  Fourier  Filters 


I.  Fourier  Analysis  of  a Periodic  Wavo 


Given  a time  dependent  periodic  function,  F(t),  it  is  possible  to  represent 

this  function  exactly  with  the  familiar  Fourier  series 

oo 

F(t)  • A^  + / (An  co 8 — t ♦ sin  ^r— - t ) (1) 

n-1 

The  sine  and  cosine  functions  of  this  series  comprise  an  orthogonal  set,  i.e., 


(a)  f0  C^i(t)  Ip  [TgCtJJn'  dt  » 0 for  all  integer  values  of  n and  n* 

PT 

(b)  L (MOL  ■ k (a  constant)  if  n«m  , 

' ' •*>  1 n l m mis  some 

-0  if  n/m  muXtiple  of  n 

2T 

(c)  /c  (f2(Oln,  (f?(0]m,  - k>  (a  constant)  if  n«  - m«  m,  ls  some 

multiple  of  n' 


where t 


if  n*/*» 


corresponds  to  cosine  terms. 


(f„(t) ] , corresponds  to  sine  terms,  and 

2 n1 

2T  - period  of  [f^(t)ln  and  [fg(01n,  for  n » n'  "1. 

Using  the  orthogonal  properties  of  equations  (2)  it  is  a simple  ratter  to 
determine  the  coefficients  ^ and  foi*  a particular  fimction  F(t)»  For 
example,  to  determine  A^  multiply’  both  siios  of  equation  (1)  by 

[f1(t)]1  - cos  t 

and  integrate  over  the  period  2T, 


2T  2T  2T 

fo  P(t)  cos  (^)t  dt  - Aq  fQ  cos  (^)t  dt  + ^ fQ  cos  (^)t  cos  (^£)t  dt 


2T 


P 2T 

.£  cos  (^t  cos  (2l)t  dt  + **•  + An  £ cos  (^t  cos  (f)t  dt 


2T  2T 

+ \ £ sin  (?)  t cos  (■^t  dt  + £ sin  (^P)t  cos  (?)t  dt  + ... 


2T 


or 


2T 


+ Bn  ^ sin  (ijriit  cos  (?)t  dt. 
2T 


.£  F(t)  cos  (^L)t  dt  - 0 + Ax  ^ cos2  (pt  dt  + 0 + ...  + 0 + 0 + 0 + ... 


+0 


According  to  equation  (2). 

2T 

^ F(t)  cos  (?)t  dt 

*1  " “2F 


• • 


^ cos  (pt  dt 


In  general. 


2T 


/ F(t)  h(t)  dt 
(a)  K'  n n 


2T  _ 

fQ  [hn(t)]2dt 


(3) 


where: 


represents  either  AR  or  Bn 
hn  (t)  represents  either  coa  (^t  or  sin  (^t. 


One  can  normalize  the  above  equation  by  requiring 

2T  i2 


^ [hn(t)r  dt  - 1 


then. 


2T 


(b)  Kn  - fQ  F(t)  hn(t)  dt  n - 0,  1,  2,  3,  ... 


(3) 


3 


where  h (t)  ■ 1 end  K • A . 
o oo 

The  normalized  Fourier  serios  can  now  be  written  as 

(c)  F(t)  - X kn  hn(t) 

n«0 

Since  the  integrals,  equations  (2b)  and  (2c),  yield  a result  differ?.'.! 
ftam  zero  only  when  the  two  functions  under  the  integrals  are  the  same,  one 
can  say  that  when  the  integral  of  (3b)  has  a value  for  a particular  hn(t), 
the  normalized  function  F(t)  must  contain  a harmnie  which  is  the  same  as 
hn(t)|  further  the  magnitude  of  Kr  is  equal  to  the  ratio  or  the  magnitude  of 
this  harmonic  wave  contained  in  F(t.)  to  the  magnitude  of  hn(t). 

II.  Finite  Fourier  Analysis 

In  practical  situations,  the  time  function  to  be  analyzed  may  not  be  a 
periodic  function.  More  often  than  not,  the  function  in  finite  with  respect 
to  time,  i.e.,  the  function  exists  only  for  a certain  duration  in  time. 

The  problem  is  to  devise  an  electrical  network  vhloh  can  analyse  finite 
functions  of  time  in  terms  of  the  cosine  and  slr.s  functions  of  a Fourier  series 
Such  a network  might  be  in  the  f orm  of  a 3 port  network  as  in  Figure  la.  Ths 

output  C would  be  proportional  to  the  integral 

T 

£ F’(t)  hn(t)  dt  (!,) 

where i 

T - half  the  period  of  h^(t) 

i 

F (t)  ■ finite  function  with  time  duration  < T. 

^ **  (*-)  exists  for  t >f,  F (t)  is  "truncated"  into  a function  F (t). 


Multiplying  C , at  tine  T by  the  appropriate  constant  would  yield  the  Fourier 
n 

coefficient  1^,  if  F (t),  a finite  function,  was  a periodic  odd  function  of 
tins,  F(t)  (equation  3b).  This  can  be  accomplished  by  making  an  odd  periodic 
extension  of  the  finite  function  f' (t)  as  shown  in  Figure  2(a). 

One  oould  form  the  3-port  network  into  a 2-port  network  by  making  h (t.) 
an  integral  generator  of  the  box  in  Figure  1(b).  The  2-port  network  can  be 
realised  as  a network  synthesis  problem  if  Cn  ia  associated  with  an  output 
function  and  F(t)  as  an  input  function.  Taking  the  ratio  of  the  Laplace  trans- 
forms of  F(t)  (the  odd  periodic  extension  of  F'(t))and  Cn,  the  ratio  of  the 
resulting  polynomials  in  s,  Q(s)  - Cn(s)  f will  allow  the  formation  of  a 

FTST 

realisable  network,  provided  Q(s)  satisfies  tho  Hurtwitz  criteria  for  a 
positive-real  rational  function  in  s^. 

Multiplication  in  the  frequency  domain,  i.c.,  F(s)  Q(s)  • Cn(a), 
corresponds  to  convolution  in  the  tim^  domain 

Cn  (x)  " °n(x)  * F(x)  " Qn(x“t)  dt  (5) 

where  t 

x ■ a particular  value  of  time  0 < x < 2T 

(^(t)  is  the  impulse  response  of  the  network  whose  function  in  the  variable  s 
is  Q(s). 

(x-t)  is  a period  of  time  for  which  (^(x-t)  is  integrated  with  F(t',  l.e., 
^(x-t)  exists  for  0 < t <x 
(^(x-t)  “ 0 for  t>x 


U3  Quillemin,  Synthesis  of  Passive  Networks,  1957,  John  Wiloy  & Sons, 
pgs.  10-36. 
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The  convolution  process  is  Illustrated  in  Figure  2.  Tho  function  0n(t), 
shown  as  a sine  wave,  is  refit  ted  about  the  vortical  time  reference.  Figure 
2(b),  and  advanced  x seconds.  The  crossod  hatched  aroa,  rhown  in  Figure  2(d), 
is  proportional  to  at  time  t“T"x.  F(t)  can  be  thought  of  as  a weighting 
function  of  the  area  under  the  curve  (^(x-t.)  as  (^(x-t)  pusses  from  x’O  to 
x-2T.  However,  wo  aro  interested  in  tho  value  at  tLno  xJT,  Figure  2(d), 
since  the  value  of  tho  integral,  equation  5,  is  proportional  to  tho  desired 
Fourier  coefficient  at  tills  instant. 

If  equation  (5)  is  compared  with  equation  (3b), 

Cn(x)  m?/r  F(t)  °n(x“t)  dt  X -T  (S) 

Kn  " •£  F(t)  hn(t)  dt  (3b) 

It  becomes  apparent,  that  over  the  range  0<t<2T,  « 1^  provided  hn(t)  - 

(^(T-t).  Wtien  h^(t)  is  a sine  function 

sin  / nTTxt  - slnfr*/  nTT\  (-t)  ) n « odd  and  ... 

'“f  f l-if — / T - TT  (&) 


• “ Cn(T)  for  0<*d  9^n3  functions. 

When  h (t)  is  a cosine  function 
n 


cos.  rvrr.  t - - cosfn*.  ntr.  (-t)] 

l ~!p  ) V ‘fl 


(7) 


• • - -(^(T)  for  odd  cosine  functions. 

Taking  into  account  the  above  results,  one  oan  substitute  Q (T-t)  for 

n 

^(t)  in  equation  (3c)  and  obtain  the  Fourier  series  of  F(t)  by  changing  the 
sign  of  the  coefficients  of  tho  odd  cosine  functions  ir.  tho  serios,  It  should 
be  noted  we  have  shown  that  the  convolution  integral  can  be  used  to  obtain  the 
Fourier  coefficients  of  an  odd  periodic  wave  at  tins  T,  where  T is  tha  half 


period  of  the  fundamental  sine  and  cosine  functions  of  the  Fourier  series. 
The  problem,  now,  is  to  construct  a network  with  ar  impulse  response 


function  Gn(t),  which  will  yield  a Fourier  coefficient  of  a given  finite 

input  function,  F'(t),  a3  an  odd  periodic  extension  F(t),  as  shown  in  Figure 

2(a),  and  analyse  only  the  first  T seconds  of  the  wave.  It  is  easy  to  show 

that  T 2T 

F'(t)  On(x-t)  dt  - ^ fQ  F(t)  0n  (x-t)  dt  (8) 

when  F(t)  is  the  odd  periodic  extension  of  F^t),  i.e.,  it  is  necessary  to 
analyse  the  finite  function  F*(t)  for  only  half  the  period  of  it3  odd 
periodic  extension  to  determine  the  Fourier  coefficients  of  the  odd  periodic 
extension  of  F'(t). 

It  should  be  noted  that  F'(t)  will  be  analyzed  as  an  odd  periodic  wave 
of  period  2T  by  the  proposed  network,  regardless  of  the  duration  of  F*(t), 
therefore  it  is  important  to  insure  that  all  significant  information  carried 
in  F'(t)  occur  within  a time  interval  such  that  F*(t)  exists  for  OSt^T. 
However,  F1 (t)  may  exist  for  an  interval  greater  than  T,  provided  the  infor- 
mation desired  lies  within  the  period  O^t^T,  such  a case  is  illustrated  in 
Figures  2(a). 

(2) 

Consider  the  two  series  connected  circuits  shown  in  Figure  3(a)  . Let 

us  assume  the  input  current  wave  Gn(t)  to  be  a sine  wave  of  period  2T.  This 
input  current  produces  a voltage  from  network  I,  and  a voltage  from 
network  2 as  shown  in  Figure  3(b)  and  3(c).  la  a unit  impulse  voltage  and 
e2  a differentiated  square  wave  voltage  of  unit  amplitude.  If  the  networks 
and  are  linear  passive  networks  a voltage  nan  be  applied 


The  method  described  here  is  from  Ouillomir.-Synthaaip  of  Passive  Networks, 
1957,  pgs.  717-726. 


r 

l 

L 

l 
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I 


i 


1 


f 


I 


across  the  sories  combination  of  and  to  produce  tho  current  wavo  Q^(t). 


Now  suppose  a voltage  o(t)  ■ o (t)  ♦ e (t-T)  la  applied  to  the  network 


combination,  Figures  3(d)  and  3(o).  One  notos  that  e(t.)  is  nothing  more  than 


a single  impulse  of  two  units  magnitude  and  that  the  output  current,  G^(t)  is 


merely  the  positive  hall  of  a sine  wave,  but  this  is  the  wave  form  nnedod  to 
determine  the  fundamental  Fourier  coefficient  A^  for  tho  finite  function  F*  (i-J • 
Therefore,  if  a network  can  be  synthosizod  with  tho  transfer  functions 

T - <?>  °„(t> 


n " If  2f  3f  ••• 


(9) 


e-j(t)  + Op  ( t) 

tho  Fourier  cooffioionto  of  F*(t)  can  bo  obtained  from  tho  networks.  Tho 

numerator  in  equation  (9)  i3  multiplied  by  2 so  that  a \ ir.it  impulse  11)  (t) 

o 

will  produce  Tn  instoad  of  2Uo(t). 

One  is  now  facod  with  synthosizing  equation  (9)  into  a roalizablo  passive 


network.  This  can  bo  accomplished  by  taking  tho  Laplace  transform  of  Tfi  to 


obtain  a ratio  of  finite  polynomials  in  tho  variablo  s.  Since  the  problem 
under  consideration  requires  Qn(t)  to  bo  tho  sine  and  cosine  terms  of  a 
Fourier  series,  it  is  best  to  express  G^(t)  an  tho  Laplace  transform  of  tho 
general  expression  for  the  Fourier  series  of  an  odd  function, 

n r n “ 


(10) 


X 2 


(A  sin-I}-^ t + P cos  -QP  t)  ■ 
n T n T 

integora  integers 


n - odd 


V * 3 \ 

n"  cild  3^ 


Tho  Laplace  transform  of  o^(t),  for  tho  purpose  of  synthesis,  can  bo  *rit.t»n 
as  the  Laplace  transform  of  tho  Four lor  series  of  a differentiated  unit  square 

w*ve,  _ 


•2<t>  • s T.  -^r— u- 

d n-  odd 


n»  odd 
integers 


<11> 


11 


(12) 


The  Laplace  transform  of  e^(t)  is  simply 

1. 


£ 


e,(t) 

Equation  (9)  can  now  be.axprossed  as 


t 


A + S 


(S)  = 


n 

“ "n 

f?  + 

(^) 

T" 

u/ 

1-  odd 

T 

(13) 


where  the  infinite  series  for  tho  square  wave  has  been  truncated  to  a finite 
number  of  terms,  J,  to  allow  synthesis  of  Tn(s). 

One  method  of  realizing  oquation  (lj)  is  to  employ  a negative  feedback 
network.  Tho  sine  functions  in  the  nomerator  of  oquation  (13)  can  be 
generated  (see  Figure  lj)  as  the  impulse  response  of  notworks  with  voltage 
transfer  functions  of  the  form 


T(S) 


(1/LC) 

S*  + (l/LC) 


(lit) 


This  transfer  function  has  tho  same  form  as  the  Laplace  transform  of  a sine  wave 
with  a frequency  of  l/LC  • 

If  a group  of  notworks  having  transfer  function  of  the  form  in  equation 
(lli),  and  time  base  outputs  which  are  odd  multiples  of  the  timo  base  2T,  are 
connected  in  parallel,  tho  sine  function  outputs  can  be  combined  in  an  adder 
such  that  the  output  voltage  will  be  a squaro  wave  (the  square  wave  in  equation 
(13)  is  a summation  of  odd  sino  wave  functions).  Tho  output  square  wave  voltage  «. — 
can  then  bo  differentiated  and  fed  back  to  the  input  as  negative  feedback,  thus 
producing  finite  sino  functions  at  tho  output  (see  Figure  ?).  Writing  the 
resultant  voltage  transfer  expression  for  tho  system  in  Figure  5» 


2 


Eout  <S>  • Eln  (S> 


* . Cs) 
T.  (a)  . -2£ 

Em  <s) 


S ) 

Y.  -Sl>(Eout(S)  £ K (S)  1 


n ■ odd 


n«  odd 


L K <s) 

n-  odd 


1 * stK  *out  (S>  f K <s>  1 

Ah.  j j 


n»  odd 


Where i 

h^(s)  ■ Laplace  transform  of  network  whose  impulse  responses  are  odd  sine  wave 
voltages* 

l{  ■ constant* 

This  is  the  form  of  equation  (13)»  except  that  we  may  also  require  cosine  functions 
in  the  numerator.  Noting  that  the  current  through  capacitor  has  the  same  wave- 
form as  the  voltage  across  the  capacitor  except  shifted  90°  in  phase,  it  is  a 
simpler  matter  to  convert  the  current  in  the  sine  function  circuits  to  cosine 
voltages.  One  can  also  obtain  the  cosine  functions  by  differentiating  the  sine 
function  outputs.  If  we  tap  off  the  sine  and  cosine  functions  from  box  1 in 
Figure  5»  they  can  be  fed  to  variable  gain  amplifiers  and  summed  to  form  a 

finite  function  F* (t)  when  the  input,  E is  an  impulse  voltage;  or  the  output 

in 

functions  from  box  1 can  be  used  to  obtain  the  Fourier  coefficients  of  the 
finite  function  F'(t)  when  Ein  is  the  function  f' ( t) . 


Ill 


III.  Description  of  the  Finite  Fourier  Filter  Unit  Proposed  by  Melpar 

Figure  6 shows  the  system  diagram  of  the  finite  Fourier  filter  unit 
proposed  by  Melpar.  The  principle  of  operation  is  the  same  as  the  system 
in  Figure  5,  section  II 5 however,  the  unit  proposod  by  Melpar  will  employ 
active  R-C  networks  instead  of  passive  R-L-C  networks,  and  the  odd  cosine 
functions  will  be  obtained  by  differentiating  the  odd  output  sine  functions. 

Ascription  of  Operation 

Odd  sine  wave  voltages  aro  obtained  from  filters  F..  to  F,  when  an  im- 

6 

pulse  is  applied  to  their  inputs.  Those  output  sine  functions  are  weighted, 
suirmod,  and  differentiated  in  block  to  form  a differentiated  square  wave. 

The  output  of  block  A,  is  fed  back  to  the  parallel  inputs  of  F,  to  F,  as 
negative  feedback.  A square  wave  generator,  block  A^,  is  used  in  conjunction 
with  the  differentiator  of  block  A^  to  provido  positive  and  negative  impulses. 
These  impulses  activate  the  filters  when  they  are  used  a3  a function  generator. 

With  the  feedback  loop  connected,  the  outputs  of  the  filters  become 
finite  functions  with  a time  duration  of  T.  Those  finito  outputs  are  tapped 
off  F^  to  F^  and  fed  to  gain  switches  1 to  6.  Each  gain  switch  contains  a 
pair  of  two  section  wafer  switches  connected  in  parallel.  The  outputs  from 
one  of  each  pair  of  switches  aro  summed  in  the  bipolar  adder,  block  A^,  to 
form  a sum  of  odd  sine  wave  functions.  The  other  outputs  from  the  gain 
switches  aro  differentiated  and  summod  in  the  bipolar  adder  to  form  a sum  of 
odd  cosine  functions.  The  odd  sine  and  cosine  functions  aro  combinod  in  block 
Ajj  to  form  a single  output.  The  unit  proposed  by  Melpar  will  also  include 
separate  outputs  for  each  sine  and  cosine  function;  thus  when  an  input  function 
is  applied  to  the  filters,  instead  of  impulses,  the  outputs  will  provide  the 
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III.  Description  of  the  Finite  Fourier  Filter  Unit  Proposed  by  Melpar 

Figure  6 shows  the  system  diagram  of  the  finite  Fourier  filter  unit 
proposed  by  Meipar.  The  principle  of  operation  is  the  same  as  the  system 
in  Figure  5,  section  II}  however,  the  unit  proposed  by  Molpar  will  employ 
active  R-C  networks  instead  of  passive  R-L-C  networks,  and  the  odd  cosine 
functions  will  be  obtained  by  differentiating  the  odd  output  sino  functions. 
Description  of  Operation 

Odd  sine  w^ve  voltages  are  obtained  from  filters  F..  to  F when  an  im- 

■i  o 

pulse  is  applied  to  their  inputs.  Those  output  sine  functions  are  weighted, 

summed,  and  differentiated  in  block  A^  to  form  a differentiated  square  wave. 

The  output  of  block  A,  is  fed  back  to  the  parallel  inputs  of  F.  to  F,  as 
• J o 

negative  feedback.  A squaro  wave  generator,  block  A^,  is  used  in  conjunction 
with  the  differentiator  of  block  A^  to  provido  positive  and  nogntive  impulses. 
These  impulses  activate  the  filters  when  they  are  used  as  a function  generator. 

With  the  feedback  loop  connected,  the  outputs  of  the  filters  become 
finite  functions  with  a time  duration  of  T.  These  finite  outputs  are  tapped 
off  F^  to  F^  and  fed  to  gain  switches  1 to  6.  F.ach  gain  switch  contains  a 
pair  of  two  section  wafer  switches  connected  in  parallel.  The  outputs  from 
one  of  each  pair  of  switches  arc  suiranod  in  the  bipolar  adder,  block  Aj^,  to 
form  a svun  of  odd  sine  wave  functions.  The  other  outputs  from  the  gain 
switches  are  differentiated  and  sunrnod  in  the  bipolar  adder  to  form  a sum  of 
odd  cosine  functions.  The  odd  sino  and  cosine  functions  are  combined  in  block 
to  form  a single  output.  The  unit  proposed  by  Molpar  will  also  include 
separate  outputs  for  each  sine  and  cosine  function;  thus  when  an  input  function 
is  applied  to  the  filters,  instead  of  impulses,  tho  outputs  will  provide  the 
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Figure  6 System  Diagram  of  Finite  Fourier  Filter  Unit 

Proposed  by  Melpar 
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individual  and  combined  sum  of  the  Fourier  coefficients  of  the  paat  T seconds 
of  the  input  function  at  time  T. 
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IV.  applications 

Prom  the  previous  sections  it  is  apparent  that  finite  Fourier  filters 
can  be  applied  to  a variety  of  analysis  and  identification  problems.  Two 
particular  applications , identification  by  modeling,  and  adaptive  matched 
filters,  have  been  studied  by  Helper  and  will  be  presented  in  the  following 
paragraphs . 

(A)  Adaptive  Matched  Filter 

It  oan  be  shown^)  that,  in  the  presence  of  white  noise,  the  maximum 
signal  to  noise  ratio  at  the  output  of  a network  occurs  when  the  network 
is  matched  to  the  input  signal.  A network  is  matched  to  a given  signal 
f(t)  if  its  impulse  response  is  f(x-t).  From  previous  considerations  it 
Was  found  that  a set  of  finite  Fourier  filters  can  generate  a desired 
finite  function  f' (t)  when  an  impulse  is  applied  to  their  inputs}  thus, 
it  is  possible  to  match  a set  of  Fourier  filters  to  a given  function  f(t) 
if  F'(t)  • f(x-t). 

Application  of  matched-filter  techniques  to  communications  and  radar 
are  well  known.  If  the  signal  structure  is  not  well  known,  in  the  sense 
that  the  transmitted  signal  has  been  distorted  in  some  unpredictable  manner, 
then  the  filters  matched  to  idealised  signals  will  no  longer  be  perfectly 
matched  to  be  actual  received  signals.  Thus,  one  anticipates  a deter ioa- 
tion  from  optimum  considerations.  Furthermore,  there  is  no  guarantee  that 
the  noise  source  will  be  white. 

The  technique  described  herein  provides  an  "adaptive'*  filter  which  tracks 
a signal's  structure  so  as  to  maintain  maximum  signal-to-noise  ratio.  This 

(3)  An  Introduction  to  Hatched  Filters,  0.  L.  Turin,  IRS  Transactions  on 

Information  Theory,  June  I960,  pp.  312. 
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is  accomplished  by  what  is  described  as  Melpar's  Adaptive  Matched  Filter 
Tracking  System. 

The  single-input,  multiple-output  linear  network,  described  in  Section 
III,  has  each  of  its  outputs  connected  to  a summing  amplifier  through  variable 
gain  input  networks.  The  response  of  this  network  can  be  made  to  approximate 
an  arbitrary  waveform  F'(t)  within  a time  interval  (0,T).  Changes  in  the 
network's  response  characteristics  are  accomplished  by  changing  the  input 
gains  on  the  summing  amplifier. 

The  value  of  the  specific  adaptive  matched  filter  described  above  stems 
from  a theorem  showing  that  each  of  the  inputs  to  the  summing  amplifier  oan 
be  independently  adjusted  for  maximum  signal-to-noise  ratio,  so  as  to  arrive 
at  a match  for  a given  input  signal.  That  is,  the  process  of  matching  the 
filter  to  a given  signal  can  be  accomplished  by  adjusting  each  input  gain 
once,  to  achieve  maximum  signal-to-noise  ratio. 

Application  of  the  adaptive  matched-filter  system  for  tracking  low-level 
signals  is  illustrated  in  Figure  7.  Signal  inputs  are  continuously  processed 
by  the  "real-time  analyzer"  network,  providing  M output  components.  These 
outputs  are  fed  to  three  parallel  circuits,  each  of  which  consists  of  a set 
of  input  networks  (controllable)  and  a summing  amplifier.  Channel  1 is  the 
matched-filter  system.  Channels  2 and  3 are  representative  experimental 
matohed-filter  systems,  where  one  of  the  input  gains  has  been  perturbed  so 
as  to  obtain  a small  increase  in  one  and  a small  decrease  in  the  other.  The 
outputs  of  channels  2 and  3 are  compared  (from  a timing  signal  generated  from 
the  output  of  channel  1)  to  establish  the  appropriate  change  in  that  input 
gain.  This  is  transferred  to  all  three  channels,  and  the  process  sequenced 
to  the  next  input. 


The  time  required  for  the  system  to  sequence  through  each  of  the  M input 
gains  is  the  same  as  the  time  required  to  receive  M inputs.  This  ~an  be 
viewed  as  the  tracking  system's  time  constant.  By  increasing  the  number  of 
experimental  channels,  it  is  possible  to  proportionately  decrease  this  time 
constant.  This  type  of  system  should  provide  improved  communications  when- 
ever signal  properties  change  in  a relatively  continuoys  fashion,  and  where 
the  signal  distortions  are  not  accurately  predictable. 

(B)  Modeling 

Modeling  is  the  technique  which  attempts  to  duplicate  a specific 
system  for  the  purpose  of  analysis.  Important  applications  of  modeling  exist 
where  it  is  impossible  or  undesirable  to  disturb  a system  which  is  in  operation. 
These  applications  are  found  in  such  areas  as  system  identification  and  experi- 
mentation for  modification  of  input  variables. 

Figure  8 shows  a block  diagram  for  modeling  a linear  system.  The  problem 
is  to  construct  a model  whose  impulse  response,  f'(t),  is  the  same  as  the  un- 
known linear  system.  Once  this  is  determined,  it  is  a simple  matter  to  write 
down  the  transfer  function  of  the  unknown  system.  The  model  consists  of  a set 
of  m Fourier  filters  whose  outputs,  Zm,  are  fed  to  a summing  amplifier  with 
variable  gain  inputs.  When  the  output  of  the  summing  amplifier  exactly  matches 
the  output  of  linear  system  the  impulse  response  of  the  unknown  linear  system 
and  the  model  will  be  the  same.  Practically,  one  can  only  hope  to  obtain  a 
good  approximation  of  the  unknown  system}  however,  the  greater  number  of  filters, 
m,  the  better  the  approximation.  It  can  be  readily  shown^  that  the  beat 
approximation  for  a given  set  of  filters  will  be  obtained  when  the  root  mean 

TUI  Melpar  Second  Triannual  Report,  W3112.01,  Machine  Intelligence  and  Adaptive 
Systems,  pp.  52-59. 
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Further,  since  the  impulse  response  is  determined  only  by  the  gain  settings 
of  the  input  amplifiers  of  the  summing  amplifier,  e can  be  used  as  a feed- 
back voltage  to  converge  the  impulse  response  of  the  model  to  f.t).  This 
implies  that  the  variable  gain  elements  in  the  model  can  be  independently 
adjusted  with  respect  to  one  another.  If  the  variable  gain  elements  remain 
independent  ovor  a specific  range,  the  modeling  system  will  be  able  to  track 
a nonstationary  linear  system  in  that  range;  providing,  of  course,  thdt the 
tracking  time  constant  of  the  model  is  negligible  compared  to  the  rate  of 
change  in  the  linear  system. 
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APreNDIX 


Demonstration  of  an  Adaptive  Matched  Filter  Breadboard 

Equipment  was  built  to  demonstrate  the  following! 
a.  The  ability  to  approximate  any  desired  finite  duration  waveform 
from  impulse  response  from  linear  networks, 
h.  The  application  of  this  concept  to  form  an  adaptive  matched  filter, 
a.  The  use  of  the  active  R-C  synthesis  technique,  recently  developed 
at  Melpar. 

The  breadboard  consists  of  two  basic  networks  which  are  essentially 
identical  in  nature.  One  of  these  functions  is  to  accomplish  the  approxi- 
mation described  in  (a)  above.  The  second  of  these  networks  corresponds  to 
the  adaptive  matched  filter  as  in  (b)  above. 

The  approximation  and  adaptive  matched-filter  concepts  have  been 
previously  reviewed.  It  corresponds,  essentially,  to  the  odd  harmonic 
Fourier  approximation.  The  approximating  terms  (impulse  responses  of  sub- 
networks) were  each  of  finite  duration  T*  In  the  breadboard,  this  time- 
duration  was  selected  as  6.25  milliseconds.  This  long  duration  was  selected 
to  illustrate  the  application  of  the  new,  active  R-C  synthesis  technique  to 
circuits  whose  L-R-C  realization  would  require  large  inductors. 

A simplified  block  diagram  of  the  breadboard  is  provided  in  Figure  A-l. 
It  illustrates  tho  design  concept  rather  than  the  design  itself.  Positive 
and  negative  impulses  (obtained  from  the  differentiation  of  a square  wave 
are  fed  to  filter  bank  A.  The  resulting  outputs  of  filter  bank  A are  the 
first,  third,  fifth,  seventh,  ninth,  and  eleventh  sine  and  cosine  waveforms 
shown  in  Figures  A-2A  through  A-2L.  It  is  noted  that  the  higher  order 
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harmonics  begin  to  deviate  from  the  Ideal,  due  to  the  absence  of  higher  order 
terms*  The  cosine  terms  appear  to  be  somewhat  more  noisy,  but  are  superior 
to  what  could  be  reasonably  expected*  Linear  combinations  of  these  response 
functions  were  used  to  generate  various  waveforms. 

The  seoond  function  of  the  equipment  is  to  provide  a filter  which  is 
matched  to  the  signals  being  generated.  Although  the  matched  filter  was 
adjusted  by  hand,  the  application  of  appropriate  circuitry  for  an  adaptive 
matohed  filter  is  obvious.  An  example  of  an  adaptive  matched  filter  system 
was  given  in  Figure  5»  section  IV.  In  Figure  A-3»  the  ninth  sine  harmonic 
was  generated  (upper  trace)  and  fed  into  a filter  matched  to  it.  The  output 
(lower  trace)  is  seen  to  be  a close  approximation  of  the  anticipated  auto- 
correlation of  the  compactly  carried  ninth  harmonic.  In  Figure  A-b,  a 
square  wave  was  generated  (upper  trace)  and  fed  into  a filter  matched  to  it 
(lower  trace)  with  the  anticipated  results. 

The  behavior  of  the  matched  filter  in  the  presence  of  noise' is  quite 
dramatic.  This  is  illustrated  in  Figures  A-5A  through  A-5D.  The  upper 
traces  correspond  to  the  signal  input  to  the  matched  filter.  For  Figures 
A-5C  and  A-5D  (upper  trace),  the  scope  sensitivity  is  2 volts/cm.  In  each 
of  these,  the  sensitivity  of  the  lower  trace  was  set  at  0.5  volt/cm.  Noise 
was  obtained  from  a "sounvister1*  white  noise  diode  and  added  to  the  generated 
signals.  Its  characteristics  are  essentially  Gaussian  from  2 cps  to  100  me. 
As  shown,  the  matched-filter  output  can  readily  detect  the  presence  of  the 
signal. 

The  above  experiment  was  repeated  with  an  asymmetric  waveform  fed  to  a 
filter  matohed  to  it  (triangular  waveform).  The  results  are  illustrated  in 
Figures  A-6A  through  A-6C.  This  was  as  anticipated. 
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Figure  4-2 
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Figure  A-3.  Ninth  Harmonic  (Sine)  - Top 

Output  of  Filter  Matched  to  the  Ninth  Harmonic  - Bottom 


Figure  A -4.  Fourier  Approximation  of  a Square  Wave  - Top 

Output  of  Filter  Matched  to  Square  Wave  - Bottom 
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SQUARE  WAVE,  1 VOLT  CM 


MATCHED  FI  LTER  OUTPUT,  0.5  VOLT  CM 

B 

SQUARE  WAVE  WITH  WHITE  NOISE,  1 VOLT  CM 
MATCHED  FILTER  OUTPUT,  0.5  VOLT  CM 

c 

SQUARE  WAVE  WITH  INCREASED  NOISE,  2 VOLTS  CM 
MATCHED  FILTER  OUTPUT,  0.5  VOLT  CM 


D 

SQUARE  WAVE  WITH  INCREASED  NOISE,  2 VOLTS  CM 
MATCHED  FILTER  OUTPUT,  0.5  VOLT  CM 


Figure  A -5 


ASYMMETRIC  WAVEFORM,  1 VOLT  CM 


MATCHED  FILTER  OUTPUT,  0 5 VOLT  CM 


ASYMMETRIC  WAVEFORM  WITH  NOISE,  1 VOLT  CM 


MATCHED  FILTER  OUTPUT,  0 5 VOLT  CM 


ASYMMETRIC  WAVEFORM  WITH  INCREASED  NOISE,  1 VOLT  CM 


MATCHED  FILTER  OUTPUT.  0.5  VOLT  CM 


Figure  a-6 


